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Abstract 

The present paper provides a general formula for the dimension of 
spline space over T-meshes using smoothing cofactor-conformality method. 
And we introduce a new notion, Diagonalizable T-mesh, over which 
the dimension formula is only associated with the topological information 
M of the T-mesh. A necessary and sufficient condition for characterization 

I I a diagonalizable T-mesh is also provided. Using this technique, we find 

that the dimension is possible instable under the condition of [1] and we 
" ^ also provide a new correction theorem. 

Keywords: Dimension, T-splines, Spline space, T-mesh, smoothing cofactor- 
conformality method 

in 

O 1 Introduction 

NURBS (Non-Uniform Rational B-Spline) is the standard for generating and 
. . representing free-form curves and surfaces, which is a basic tool for using in 

CAD [2] and is also a desirable tool for iso-geometric analysis ^Sj . A well-known 
and significant disadvantage of NURBS is that it is based on a tensor product 
}J] structure with a global knot insertion operation. It is desirable to generalize 

C3 NURBS to spline space which can hander hanging nodes. 

Spline spaces over T-meshes S{di,d2,a, f3,T) were first introduced in 
which is a bi-degree (^1,^2) piecewise polynomial spline space over T-mesh T 
with smoothness order a and /3 in two directions. Spline space over T-meshes 
have been applied in fitting [5], stitching simplification [7j, isogeometric 
analysis [5], [H], solving elliptic equations [TD] and spline space over triangula- 
tions with hanging nodes [11'. Spline over T-mesh has several advantages. For 
example, it has a simple local refinement which will never introduce additional 
refinement according to the definition. And it is a polynomial in each face which 
has simple and efficient integration for numerical analysis. 
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NURBS with local refinement is discovered before spline over T-mesh, called 
"T-spline", which are defined on a T-mesh by certain collections of B-splines 
functions defined on the mesh [12] , [E] . T-splines have been proved to be a pow- 
erful free-form geometric shape technology that solves most of the limitations 
inherent in NURBS. T-splines have several advantages over NURBS such as lo- 
cal refinement [TB] and watertightness [TJj , and they are forward and backward 
compatible with NURBS. These capabilities make T-splines attractive both for 
CAD and also desirable for iso-geometric analysis [15]. However, it is very diffi- 
cult to unravel the mystery of T-spline spaces which have important implications 
in establishing approximation, stability, and error estimates |16j . And till now, 
only a sub-class of T-spline space, analysis-suitable T-splines space [T7], [T5] . 
[T^ has been discovered using the technology from spline space over T-meshes. 

Thus, it is very important to understand the spline space of piecewise poly- 
nomials of a given smoothness on a T-mesh, which foundation but non-trivial 
step is to calculate the dimension of the space. Till now, many different meth- 
ods have been applied to tackle these issues, such as B-net g|, minimal deter- 
mining set [50], smoothing cofactor-conformality method [I] and homological 
technique [3T]. B-net and minimal determining set methods are suitable for 
spline space with reduced regularity, i.e., the degrees for the polynomial is larger 
enough than the smoothness order. Smooth cofactor-conformality method [22] , 
[23] is a power tool for calculating the dimension of spline space in multi-variate 
splines theory. It transfers the smoothness conditions into algebraic forms and 
calculate the dimension using linear algebraic tools. Homological technique is 
another power tool for calculating the dimension using the similar idea except 
regarding the smoothness conditions as the kernel of a certain linear maps. The 
present paper is focusing on smoothing cofactor-conformality method. 

All the existing dimension results are forcing the T-mesh to be special nesting 
structure, such as hierarchical, regular T-subdivision and no cycles [20]. With 
such structure, one can calculate the dimension level by level. Our approach is 
based on a different point of view. As smoothing cofactor-conformality method 
can convert the smoothness conditions into algebraic forms, so we directly focus 
on the algebraic forms and study the condition, under which the dimension of the 
linear system doesn't associated with knot values. And then wc use the condition 
to find a new notion, diagonalizable T-meshes, is the corresponding T-meshes. 
We believe that this new notion is the key condition to compute the dimension 
using smoothing cofactor-conformality method. The main contribution of the 
present paper includes, 

• We provide a general formula for the dimension of spline space over any 
regular T-meshes without holes using smoothing cofactor-conformality 
method. 

• We provide a new notion, diagonalizable T-mesh, over which the dimen- 
sion formula is only associated with the topological information of the 
T-mesh. We also provide a necessary and sufficient condition for charac- 
terization diagonalizable T-meshes; 
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• We discover that the dimension is possible instable under the condition 
of [Tj and we also provide a new correction theorem. 

The remaining paper is structured as follows. Pertinent background on spline 
space over T-mesh is reviewed in Section 2. In Section 3, we review how to use 
smoothing cofactor-conformality method to analysis the dimension of spline 
space over T-meshes. In section 4, we provide a condition under which we can 
calculate the dimension of spline space over T-mesh without considering the 
knot values. In section 5, we first show that the dimension result in [1] is not 
correct and we provide a new correction based on this new technology. The last 
section is conclusion and future work. 

2 T-mesh and spline over T-mesh 

In this section, we briefly review the notion of spline spaces over T-meshes and 
the related dimension results of the corresponding spline space. 

2.1 T-mesh 

A T-mesh T is a collection of axis-aligned rectangles Fi such that the interior 
of the domain is UF^, and the distinct rectangles F.^ and Fj can only intersect 
at points on their edges. The rectangles Fi are also called the face or cell of 
the T-mesh. The vertices of the rectangles are called the nodes or vertices for 
a T-mesh. The line segment connecting two adjacent vertices on a grid line 
is called an edge of the T-mesh. T-meshes include tensor-product meshes as 
a special case. However, in contrast to tensor-product meshes, T-meshes are 
allowed to have T-junctions, or T-nodes, which are vertices of one rectangle 
that lies in the interior of an edge of another rectangle. The domain J7 need 
not be rectangular, which may have holes, concave corners. For example the 
T-mesh in Figure [l] the grey region is a hole and vertices V38 and Vss are both 
concave corners. In the present paper, we require the T-meshes to be regular 
and without holes. Here regular means that the set of all rectangles for a 
T-mesh containing a vertex has a connected interior [21] . 

The vertices, edges can be divided into two parts. If a vertex is on the 
boundary grid line of the T-mesh, then is called a boundary vertex. Otherwise, 
it is called an interior vertex. If both vertices of an edge are boundary vertices, 
then it is called a boundary edge; otherwise it is called an interior edge. An 
l-edge is a line segment which consists of several interior edges. It is the longest 
possible line segment, which interior edges are connected and two end points 
being T-junctions or boundary vertices. 1-edges have three different classes. If 
the two end vertices of a l-edge are interior vertices, then the l-edge is called 
interior l-edge. If two end vertices of a l-edge are both boundary vertices, then 
the l-edge is called a cross-cut. Otherwise, if one end vertex is boundary vertex 
and the other is interior vertex, then the l-edge is called a ray. A mono-vertex is 
the intersection vertex of an interior l-edge and a cross-cut or a ray and a free- 
vertex is the intersection between cross-cuts and rays. For example, in Figure[l] 
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Figure 1: A T-mesh. 



vertices V53, V44 and V^q are interior vertices, and V15, Vig and Vgs are boundary 
vertices. The 1-edge V15V95 is a cross-cut, while V16V46 is a ray, and V44V84 and 
V57V58 are interior l-edges. 

For later use we introduce some notations for a T-mesh as shown in Table [T] 



2.2 Spline space over T-mesh 

Given a T-mesh T € M^, let denote all the cells in T and fl the region 
occupied by all the cells in T. The bi-degree ((ii,c?2) polynomial sphne space 
over T-mesh T with smoothness order a and /? is defined as 

5(di,d2,a,/3,r) := {f{x,y) G e Pdid.,V0 G J"}, 

where P^ida is the space of all the polynomials with bi-degree (^1,(^2), and 
C"'''(r2) is the space consisting of all the bivariate functions which are contin- 
uous in n with order a along x direction and with order /3 along y direction. 
It is obvious that iS(di, (^2, a, /3, T) is a linear space, which is called the spline 
space over the given T-mesh T. 

Until now, several articles have been studied to analysis the dimension of 
the spline space over some special families of T-meshes. 

• Reduced regularity: 

In 2006, [4j studied the dimension of the spline space under the constrains 
that the order of the smoothness is less than half of the degree of the spline 
functions. According to Theorem 4.2 in |;4 , it follows that if c?i > 2q! -I- 1 
and d2 > 2/3 + 1, 

dimS{di,d2,a, /3,T) =F{di + l)(d2 + 1) - -B^K + + 1)- 
Ey{d2 + l){a + 1) + V{a + 1){I3 + 1), 
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Table 1: Notations for a T-mesh 



F 


number of faces in T 




number of norizonal interior edges m / 


rpv 


number oi vertical interior edges m / 


V 


number of interior vertices in T 


fjh 


number of horizonal cross-cuts in T 


QV 


number of vertical cross-cuts in T 




number of horizonal interior 1-edges in 7" 




number of vertical interior 1-edges in T 


rie 


number of interior 1-edges in T (T^ -f T") 


v+ 


number of free-vertices in T 


N''- 


minimal integer larger or equal to ^^'-^^ 




minimal integer larger or equal to 


ric 


(di-a)(d2-/3)(T/- 


rir 


(dl + l)(d2 - /3)T'' + (d2 + l)(rfi - ol)T'" 



where F, E^, Ey, and V are defined in Table [T] ^24j also proved this result 
using minimal determining set method. And later, |25) analysis a special 
T-spline with reducing regularity using the dimension in 0]. 

• Enough mono-vertices 

In 2006, [1 calculated the dimension of spline space over a T-mesh if each 
interior l-edges have enough mono-vertices. In the T-mesh, if the interior 
of each horizontal interior l-edge has at least N''' — 2 mono-vertices and 
the interior of each vertical interior edge segment has at least iV" — 2 
mono-vertices, then the dimension of spline space over the T-mesh is, 

dim5(di,d2,a,/3,r) =(rfi + l)(rf2 + 1) + [Ch - n){di + l){d2 - /?) + 
(a - Ty){d2 + - a) + V{di - a){d2 - 

here Ch, Cy, T^, Ty and V are defined in Table [l] However, we will show 
that the condition in this paper is not right. 

• Instability: 

In 2011, [26] [27j discovered that the dimension of the associated spline 
space is instability over some particular T-meshes, i.e, the dimension is 
not only associated with the topological information of the T-mesh but 
also associated with the geometric information of the T-mesh. 

• Analysis-suitable T-splines: 

In 2011, [T71 [TH] provided a mildly restricted subset of T-splines, which 
optimized to meet the needs of both design and analysis. And JJ^j compute 
the dimension of the spline space S{d,d,d — l,d — 1, T) if the T-mesh T 
is an extended T-mesh of an analysis-suitable T-mesh. 

• Regular T-subdivision: 

|21j studied the dimension for spline space S{di,d2,a, I3,T) when the 
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T-mesh is a regular T-subdivision by exploiting homological techniques, 
which is a special case of the present paper. 

• Special hierarchical T-mesh: 

|28j provided the dimension for spline space S{d, d, d ~ I, d — 1,T) over a 
special hierarchical T-niesh using homological algebra technique. 



3 Smoothing cofactor-conformality method 

In this section, we will review one of the main methods, smooth cofactor- 
conformality method introduced in [22 and 23^, for computing the dimension 
of spline space over T-meshes. 

In the theory of multi-variate splines, in order to calculate the dimension of a 
spline space, one first needs to transfer the smoothness conditions into algebraic 
forms. 

Referring to Figure^ for any interior vertex Vij — {xi, i/j), suppose the four 
surrounding bi-degree (3i, 1^2) polynomial patches are Pi j{x, y),k = 0, 1, . . . , 3 
respectively (if the vertex Vi,j is a T-junction, then some of the polynomial 
patches are identical). For example for the left T-junction in Figure[2]b, patches 
Pij(a;,y), Plj{x,y) are identical. 



Pij(x,y) 




Pij(x.y) 


Hi 
















Pi!j(x,y) 











\i 

















Figure 2: smoothing cofactors around a vertex. 

As p^j{x,y) and plj{x,y) are C" continuity, so there exists a bi-degree 
(di — a — 1, rf2) polynomial Xf j{y) such that 

Phix,y) ~plj{x,y) = Xljix,y){x~x^)°'+\ (1) 

Here j{x, y) is called edge cofactor for the common edge of patches p^ j{x, y) 
and pI j{x, y). If two patches are identical, then the edge cofactor is zero. 

Similarly, there also exist bi-degree (di — a — 1,(^2) polynomial X}j{x,y), 
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bi-degree (di, ^2 — /? — 1) polynomials IJ^\ j{x, y) and /if ^(x, y), such that 

pI]{x, y) - y) = /^^.(a;, y)(y - 

Plj{x,y)~plj{x,y) = -XlJ{x,y){x~x^)"+\ 



(2) 
(3) 
(4) 



Sum with all these equations, we have 

iXl{x,y) - Xlix,ymx - x^^' = if^lix,y) - ^^Ux,y))iy - y,f+\ (5) 

Since (x — Xi)""*"^ and (y — yj)^~^^ are prime to each other, so there exist 
bi-degree (di — a — 1, (i2 — /? — 1) polynomial di,j{x, y), such that, 

•^ij(2;,2/)-A-j(a;,2/) = d,j{x,y){y-yjf+^, fil^{x,y)~filj{x,y) = d,,j(a;, 

Let 

di-a-l (i2-;3-l 

d,.,(x,2/)= ^ ^ df]{x-x,ny-y,y, 
p=0 q=a 

We call these df'j vertex cofactor. Denote to be a vector contains all 
coefficients d^'J , i.e., 



a+l 



J / ,0,0 ,0.1 



,di-Q-l,(i2-/3-2 idi-a-1, (12-/3-1 
" 1 ' 7 



)• 



For boundary vertices, there is a little different. Since in this case, we only 
have parts of the four equations such as ([I]) , ([2]) , (|3| , Q . So we don't need to 
assign the bi-degree (di — a — 1, c?2 — /3~ 1) polynomial dij{x,y) for the boundary 
vertex. Instead, we will assign the edge cofactor for the corresponding edge. For 
example, in FigurejS^, we need two edge cofactors X] j and fi] ^ and for the figure 
b, we need only one edge cofactors iijj. 



Pij(x,y) 



A,!, P,^j(x,y) P,yx,y) 



Pi,j(x,y) 



U 



Pij(x,y) 




Figure 3: smoothing cofactors around a boundary vertex. The grey regions are 
outside of the T-mesh. 
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(a) 


a. 
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Figure 4: smoothing cofactors along a horizonal edge segment. 



The vertex cofactors for the interior vertices are not totaUy free, there are 
other constrains for the continuity condition along each l-edge in the T-mesh. 

We first consider a horizontal interior l-edge referring to Figure|4]with fc + 1 
vertices i^i^.j , t = 0, 1, . . . , fc. According to equation ([6|, we have 

fJ'Hji^^y) - f^itji^^y) ^ K,j{x,y){x-x^^)°'+'^, (6) 

and 

Sum all these equations, we have the following equation, 

k 

^d,„,(x,y)(x-x,,r+i = o. (7) 

t=o 

Similarly, the constrains for a vertical interior l-edge is 

I 

Y,d,.,,{x,y){y-y,y+^^0. (8) 
t=o 

The above equations are called edge confer mality conditions. 

Lemma 3.1. Assume each Xi^ and yj^ are distinct, then the dimensions of 
solution space of equation ([7| and ([8| are (^2 — (3){l{di — a) — di — 1)+ and 
(c?i — a)(fc(c?2 — P) ~ d2 ~ 1)+ respectively. 

Proof. See [1 for more details. □ 

Remark 3.2. // the number of the vertices in a horizontal l-edge is less than 
, or the number of the vertices in a vertical l-edge is less than N"" , then the l- 
edge will not contribute the dimension of the .spline space, i.e., we can delete the 
l-edge without altering the spline space. Thus, we called .such l-edge vanished 
l-edge. In the following, we assume that all l-edges are not vanished l-edges. 
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Now we consider a horizontal ray with r + 1 vertices Vit,j,t — 0, 1, . . . , A:. 
Without loss generalization, we assume i^^-^j is a boundary vertex. According 
to equation we have 

^^i,3i^^y) ~ f^lji'^^y) ^ '^H,jix,y){x - XiX^\ (9) 

and 

fJ-l-i.ji^^y) = fJ'l,j{x,y). 
Sum all these equations, we have the following equation, 
fc 

J2d^,A^,y){x~x,X+'=fil^^ix,y). (10) 
t=i 

The main different of the constraint from ([7| and Q is that these is one edge 
cofactor fif_^ j{x, y) in the constraints. For any di^_j{x, y) assigned for the interior 

vertices, we can assign iif_^ j{x,y) as '^it.jC^^i ~ ^u)"'^^ ^o satisfy the 

constraint. 

For a horizontal cross-cut in the T-mesh, since it has two boundary vertices, 
we can conclude that it has (di + l){d2 — (3) degrees of freedoms. Similarly, a 
vertical cross-cut has (^2 -I- l)(rfi — a) degrees of freedoms. 

The linear systems as ([t]) or ([s]) , associated with each interior l-edges can be 
assembled into a global system as A^x — 0. Here M is a Ur x ric matrix, which 
is called conformality conditions matrix. And x is a column vector whose 
elements are all the vertex cofactors for the interior vertices in the T-mesh. 

And according to the above anlaysis, we can get the dimension for spline 
space over any general T-mesh without holes according to smoothing cofactor- 
conformality method which is stated in the following theorem. 

Theorem 3.3. Given a T-mesh T which has no vanished l-edges and holes, let 
matrix M. he the conformality conditions matrix, then the dimension of spline 
space over the T-mesh is, 

dim5(di, d2, a, /3, T) + l)(rf2 + 1) + C'\di + l)(rf2 - /3) + 

^(da + l){di - a) + V+{dx - a){d2 -(3) + dimTW, 

4 Diagonalizable 

Theorem |3.3| indicates that the main difficult to compute the dimension of spline 
space over T-mesh is to calculate the rank of conformality condition matrix Ai. 
It is obvious that the structure of the matrix is associated with the order of 



9 



edge conformality conditions and it is also associated with the order of vertices 
cofactors. Most of existing method study the dimension of the matrix by forcing 
the T-mesh to be nest structure, such as T-subdivision, hierarchical. However, 
nesting structure should be fine for application, but it is hiding some essential 
properties for T-mesh. What we wish to answer is that under what condition 
can we compute the dimension regardless the knot intervals for a given T-mesh 
and why such condition is important. 

Definition 4.1. Given a T-mesh T , suppose we order the all the interior l- 
edges as ei.,j — 1,2, .. . ,ne, then we can compute the new-vertex-vector v"^ . 
Here v\ is the number of vertices on l-edge e^^ and Vj is the number of vertices 
on l-edge e^. hut not on l-edges e^j., /c < j. 

Definition 4.2. A T-mesh is called Diagonalizable if there exists an order 
of l-edges Ci.^j — 1,2, ...,rte such that the new-vertex-vector satisfies that 
v^j > A'^^ if Ci- is horizonal and v'^ > N'" if Ci- is vertical. 

Lemma 4.3. // a T-mesh is diagonalizable, then the matrix A4 has full column 
rank regardless the knot intervals. 

Proof. Since we assume that the T-mesh has no vanished l-edges, the matrix 
A4 has more columns than rows, i.e., Uc < n,-. After arranging the order of 
edge conformality conditions and the order of vertex cofactors, an appropriate 
partition of the linear system of constraints is 

= (11) 

where M.i is a n^. x matrix and A^2 is a x (ric — n^) matrix, is a 
vector of the first vertex cofactors, and X2 is a vector of the remaining vertex 
cofactors. 

Since the T-mcsh is diagonalizable, so there exist the order for the interior 
l-edges satisfy the condition stated in definition |4.2| Without loss of generaliza- 
tion, we assume the order for the l-edges are e^, i = 1, 2, . . . , rig- Then we arrange 
the order of the edge conformality conditions corresponding to l-edge from 
to ei. The order of the vertex cofactors can be placed in the following fashion. 
For l-edge e„^ , if it is a horizonal l-edge, then there exist vertices which are 
not appeared in the other l-edges. Each vertex corresponds [di — a) x (^2 — /3) 
cofactors, so these vertices correspond {di — a) x (^2— /?)-/V'' > {di + l){d2— P) 
cofactors. Select any (di \){d2 — /?) cofactors and put in the beginning of x|". 
If the l-edge is a vertical l-edge, we will select {d2 + l)(rfi — o) cofactors and 
put them in the beginning of x^. These process can be applied for the remain- 
ing l-edges. And then the matrix A^i is in upper block triangular form and 
according to Lemma [3^ each diagonal block (^2 — I3){di -f 1) x (^2 — /3)(c?i -I- 1) 
or (c?i — a){d2 -I- 1) x {di — a){d2 -t- 1) matrix is full rank, thus matrix A4i is 
obviously of full rank. □ 



[Mi\M2] 



Xl 
X2 
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Theorem 4.4. Suppose T-mesh T is diagonalizable and has no vanished l-edges 
and holes, then the dimension of spline space over the T-mesh is, 

dim5(di, d2, a, /3, T) + l)(d2 + 1) + {C' - T''){di + l){d2 - 

(C - r'')(d2 + l)idi -a) + - a)(d2 - /?). 

Proof. As T-mesh T is diagonalizable, so 

dimTW = = {di-a){d2-l3){V-V+)-{di+l){d2-P)E'''-{d2+l){di-a)E'' . 



According to theorem 3.3, we have the dimension dimiS(c?i, ^2, a, /3, T) is 



dim ={di + l){d2 + 1) + + l)(d2 -I3) + C^d2 + l)(di -a) + V+{d, - a){d2 

(di - a){d2 - I3){V - V+) - (di + l)(d2 - I3)E'' - {d2 + l)(di - a)E'' 
= idi + l)(rf2 + 1) + {C' - T'')(di + l)(d2 - /?) + - T''){d2 + l)(di - a)+ 
y(di-a)(d2-/3). 

□ 

Now we provide several examples for dimension of spline space over the 
following T-meshes. 











ei 









(a) a. (b) b. 

Figure 5: Two example T-meshes for sphne space 5(3, 3, 2, 2, T). 



Example 4.1. The first examples are spline space 5(3,3,2,2, T) over the two 
T-meshes in Figure^^ 

The first T-mesh has some concave corners and one interior l-edge which is a 
vanished l-edge since it has only two vertices. And the T-mesh has two cross-cut. 
So according to theorem \ 3.3[ the dimension of the spline space 5(3,3,2,2,7") 
over the first T-mesh is 16 -f 2 x 4 = 24. 

The second T-mesh is much more complex. It has four interior l-edges 
ei,i = 1, ... ,4, four cross-cuts and five rays. If we arrange l-edges as Ci, 62, 
63, 64, then the new-vertex-vector corresponds this order is (5,5,4,3). But we 
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can arrange the l-edges as order ei, 64, e^, 62, then the new-vertex-vector corre- 
sponds this order is (5, 4, 4, 4), which means the T-mesh is diagonalizable. In the 
T-mesh, we have 4 cross-cuts, 2 horizonal and 2 vertical interior l-edges and 31 
interior vertices, so according to theorem \4-4\ ihe dimension of the spline space 
5(3, 3, 2, 2, T) over the first T-mesh is 16 + 4(4 - 4) + 31 = 47. We should men- 
tion that non of existing method can calculate the dimension for this example 
since the T-mesh is not T-suhdivision or hierarchical. 



n — r 



— n 



n — r 



=J=« 



Vss 



£3 



ei 



i — « 



n — r 



i — » 



Figure 6: Non-diagonalizable vs. diagonalizable for spline space 5(3,3, 1, 1,T). 



Example 4.2. The second example is associated with spline space 5(3, 3, 1, 1, T) 
over the two T-meshes in Figure^ 

The first T-mesh has four interior l-edges. And we can see that it is not 
diagonalizable since the new-vertex-vector could he (3,2,2,1) or (3,3,1,1) for 
different orders. Thus, we cannot complete the dimension using the current 
method. Actually, according to our knowledge, no existing method can compute 
the dimension of such spline space. 

The second T-mesh also has four interior l-edges, but according to our knowl- 
edge, no existing method can compute the dimension of such spline space. But 
if we arrange the order of the l-edges to he ei, e^, 63 and 64, we can see that 
the new-vertex-vector is (3,2,2,2), i.e., the T-mesh is diagonalizable. Since 
the T-mesh has 8 cross-cut, 4 interior l-edges, and 27 interior vertices. So 
the dimension of the spline space 5(3, 3,1,1, T) over the second T-mesh is 
16 + 4 X 2 X (8 - 4) + 27 X 4 = 156. 

Example 4.3. The third example is associated with spline space 5(3, 3, 2, 2, T) 
over the two T-meshes in Figure^ Both T-meshes have four interior l-edges 
and we can check that the T-meshes are not diagonalizble. For the first T-mesh, 
the dimension is 65 and the dimension for the second T-mesh is instable, i.e., 
the dimension is associated with the value of the knots. 

We can see that the second T-mesh can be constructed by moving two red 
l-edges of the first T-mesh. Although the two T-meshes have different structure, 
hut if we look at the structure of the conformality conditions matrix, we can see 
that the structure of two matrixes are identical, except the two red l-edges using 
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Figure 7: The dimension of bi-cubic spline space over the T-mesh is instable, 
which is generated from the left T-mesh by moving two red l-edges. 



the knots of S2 to Sg instead of Si to S5. So diagonalizhle is a concept for the 
structure of the conformality conditions matrix not T-mesh itself. 

This example also tells us that if a T-mesh is not diagonalizahle, then we 
have to consider the knot values in order to analysis the dimension of the spline 
space using smoothing co factor- conformality method. 

Remark 4.5. A similar result has been abstained in '21]l, which provides the 
dimension for a special T-mesh, called regular T- subdivision. The main differ- 
ence between these two method is that the diagonalizahle T-meshes don't need 
to be nested structure. Regular T-subdivision is a special case of diagonalizahle 
T-mesh. 



4.1 Characterization 

In this section, we will provide a necessary and sufficient condition for charac- 
terization a diagonalizable T-mesh. 

Lemma 4.6. A a necessary and sufficient condition for a T-mesh to be diag- 
onalizable is for any interior l-edges set S, there at least exists one horizonal 
l-edge such that the number of vertices on this l-edge but not on the other l-edge 
in S is at least , or there at least exists one vertical l-edge such that the 
number of vertices on this l-edge but not on the other l-edge in § is at least N"" . 

Proof. First, we prove the condition is necessary using reduction to absurdity. 
IF the T-mesh is diagonalizable, but there exists a set of l-edges {e^ ^ , e^^ , . . . , } 
such that any horizonal l-edges in the set at most have vertices which are 
not on the other l-edges in the set and any vertical l-edges at most have TV" 
vertices which are not on the other l-edges in the set. And since the T-mesh 
is diagonalizable, so without loss of generalization, we assume when the order 
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of the interior 1-edges is 61,62, .. . ,6„_^, in which order the 1-edges satisfy the 
condition of diagonahzable. Let k to the maximal index for all = 1, . . . ,s. 
Now, we consider 1-edge 6fc, since it has at most A^'' — 1 or — 1 vertices which 
are not on the other l-edges in the set, so it also has at most A^'' — 1 or iV" — 1 
vertices are not on the 1-edges for 61, ... , Ck-i since S C {61, . . . , efe_i, 6^}, which 
violates the assumption of diagonahzable. Thus, the condition is necessary. 

Now we prove the condition is sufficient. For the set of 1-edges ei,i = 
1, . . . , ne, according to the assumption, there exist one 1-edge which has enough 
vertices on the 1-edge but not on the others. Without loss of generalization, we 
assume it is 61. Suppose we have ordered the 1-edges as 61, 62, . . . , e^- satisfy the 
diagonahzable condition, then for set {cj+i, . . . ,6„^}, according to the assump- 
tion, there exist one 1-edge which has enough vertices on the 1-edge but not on 
the others. Without loss of generalization, we assume it is ej+i. With this pro- 
cess, we can order the 1-edges such that it satisfy the diagonahzable condition, 
which completes the proof. □ 

5 Correction for [T] 

In this section, we will show that the dimension result in [1 is not right. Pre- 
cisely, the dimension under the condition of [T] is possible instable. And we 
also provide a new modified theorem using the technology builded in the last 
section. 

5.1 A instable example under the condition of [T] 
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Figure 8: Counterexample for paper 

Consider the spline space 5(3, 3, 2, 2,T) over the T-mesh illustrated in Fig- 
ure [H There are four interior 1-edges in the T-mesh and each one have two 
mocro- vertices in the interior, which satisfy the condition in [I]. However, we 
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will analysis the dimension in the following and show that the dimension of the 
spline space over the T-mesh is instable. 

Actually, we arrange the interior l-edges as the order of ei, 62, 63 and 64. 
And we arrange the order of the vertices as 112.2, V3,2i ^4,2, ^'8,2, V7.2, ^7.3^ vt^4, 
V7,8, i'7,7, Vqj, vij, V2,7^ «2,6, ^2,5, «2,i and W2,2, then we can get the sparse 
matrix A4 which has the following form, 
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11 Si — i and tj — j, we can verify that the rank of the matrix is 15, thus the 
dimension of the spline space is 49. If we perturb one of the knots a little bit, 
such as S3 = 3.0 + e, where e is an arbitrary small value, then the dimension is 
48. In other words, the dimension is instable. 

5.2 Correction theorem 

In this section, we will provide a new theorem to correct the theorem in J] using 
the method in the last section. 

Theorem 5.1. Given a regular T-mesh T without holes, if each horizontal l- 
edges has at least N'^ — 1 mono-vertices and each vertical l-edge has at least 
N'" — 1 mono-vertices except two end vertices, then the dimension of the spline 
space defined on T is 

dim5(di, d2, a, /3, T) =(di + l)(d2 + 1) + {C" - T''){di + l){d2 - /3)+ 
{C-" - T''){d2 + l)(di - a) + - a){d2 ~ /?). 

Proof. We will prove that under the condition, the T-mesh is diagonalizable 
using reduction to absurdity. Actually, if the T-mesh is not diagonalizable. 
Then for any set of l-edges, it is bounded. Suppose is the most bottom 
horizonal l-edges in the set (if there are more than one l-edge in the set, we 
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6 

Figure 9: The new correction theorem for [1,. 

will pick the leftmost one) , see Figure [9] as an illustration. According to the 
assumption, has at least N'^ — 1 mono-vertices (black rectangle vertices in 
the figure), if one of the two end vertices is not on the other l-edges of in the 
set, then at least has N'^ vertices which are not on the other l-edges in the 
set. According to Lemma [4. 6[ the T-mesh is diagonalizable which violates the 
assumption. Thus, there exists a vertical 1-edge, Cj, which contains one of the 
end vertices of e^. Without loss of generalization, we assume the right end 
vertex is on the 1-edge Cj in the set. Now, we consider 1-edge e^, the bottom 
end vertex of the 1-edge cannot lie on the other l-edges in the set because it is 
on the bottom of which is the bottommost l-edges in the set. So Cj at least 
has A^" vertices which are not on the other l-edges in the set. According to 
Lemma |3.1[ the T-mesh is diagonalizable which violates the assumption. Using 
Theorem |4.4[ we prove the theorem. □ 

6 Conclusion and Future work 

In the present paper, we introduce a class of T-meshes, diagonalizable T-meshes, 
over which the dimension of spline spaces is stable. We also provide a necessary 
and sufficient condition to characterize this class of T-meshes. The dimension 
result in the present paper can cover all the existing dimension results as special 
cases. 

The paper leaves several open problems for further research. As we have 
provided the dimension of the spline space, so there are many problems which 
need to be solved, such as construction of a set of basis functions with good 
properties, geometric operations and properties of the splines over T-meshes, 
etc. We will explore these problems in detail in future papers. It is also an 
important and interesting question to find out other general with fix dimension 
T-meshes. 
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